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In this work we investigate 1+2 reactions within the framework of the hyperspherical adiabatic 
expansion method. To this aim two integral relations, derived from the Kohn variational principle, 
are used. A detailed derivation of these relations is shown. The expressions derived are general, 
not restricted to relative s partial waves, and with applicability in multichannel reactions. The 
convergence of the /C-matrix in terms of the adiabatic potentials is investigated. Together with a 
simple model case used as a test for the method, we show results for the collision of a 4 He atom on 
a 4 He2 dimer (only the elastic channel open), and for collisions involving a 6 Li and two 4 He atoms 
(two channels open). 
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I. INTRODUCTION 

Calculation of phase shifts (or the /C-matrix) for a 
given reaction is often complicated by the necessity of 
knowing the wave function of the full system at large 
distances. Extraction of the phase shifts can be in prin- 
ciple achieved by comparison of the large distance part of 
the wave function with its known analytic asymptotic ex- 
pression. For processes involving only two particles (1+1 
collisions) this procedure can be easily implemented, and 
therefore the phase shifts can be computed. However, the 
more particles involved in the reaction the more difficult 
the calculation of an accurate wave function at large dis- 
tances, or at least the more expensive from the computa- 
tional point of view. Therefore, when increasing the num- 
ber of particles the extraction of the phase shifts becomes 
progressively more and more complicated. In nuclear 
physics, collisions involving three and four nucleons have 
been extensively studied solving the Faddeev (A = 3) and 
Faddeev-Yakubovsky (A = 4) equations [l|, [|| , and the 
Hyperspherical Harmonic (HH) expansion in conjunc- 
tion with the Kohn Variational Principle (KVP) [1, 0]- 
These methods have been tested through different bench- 
marks [H, 0] . When the interaction between the particles 
presents a hard core, as in the case of the atom-atom in- 
teraction, a direct application of these techniques could 
be problematic. The Faddeev equations has been mod- 
ified to deal with a hard core repulsion [tJ and, in the 
case of the HH expansion, a correlation factor has been 
included Q. In addition the Hyperspherical Adiabatic 
(HA) expansion method has proven to be a very efficient 
tool 0. 

In the case of atom-atom interactions, the HA expan- 
sion shows a particularly fast range of convergence in the 



description of bound states, as has been shown for ex- 
ample in Ref. fiol ] for the description of rare gas trimers. 
In the past years there was a systematic use of the HA 
expansion in the description of three-atom systems in 
the ultracold regime (see for example Refs.fTll. [l2| and 
references therein). These applications rise the question 
about the convergence properties of the HA method for 
scattering states, in particular in the description of a 1 + 2 
collision. In principle the HA expansion could be ap- 
plied to describe such a process since it leads to a clean 
distinction between all the possible incoming and out- 
going channels. However, as was recently showed, the 
convergence of the expansion slows down significantly in 
applications directed to describe low energy scattering 
states 13]. This problem appears at the moment of ap- 
plying the proper boundary conditions to the hyperradial 
functions. In fact, in the HA expansion, the hyperradial 
functions are obtained solving an infinite system of equa- 
tions in the hyperradial variable p and the convergence 
of the expansion is studied by increasing the number of 
equations considered after truncation of the system. For 
describing a 1 + 2 collision, the hyperradial functions are 
obtained requiring an hyperradial plane wave behavior 
as p — > oo. However, in such a process, the plane wave 
behavior results in the relative distance between the in- 
cident particle and the center of mass of the two-body 
bound system. The equivalence between both descrip- 
tions happens at p m oo or, in other words, by including 
a very large number of hyperradial functions in the solu- 
tions. This is the cause of the extremely slow observed 
convergence. 

In Ref. [3] the authors introduced a general method 
to compute the phase shift from two integral relations 
that involve only the internal part of the wave function. 
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This method is a generalization to more than two par- 
ticles of the integral relations given in [l^, [l6| and it is 
derived from the KVP. In the case of the HA expansion, 
in Ref. 14] was shown that for a 1+2 reactions, the use of 
the integral relations allows to determine the phase-shift 
with a pattern of convergence similar to a bound state 
calculation. Therefore, thanks to the integral relations, 
the hyperspherical adiabatic expansion method appears 
as a powerful tool also to describe scattering processes. 

The purpose of this work is to show in details the use 
of the integral relations in conjunction with the HA ex- 
pansion method to describe scattering states. In Ref. [l4| 
the particular case of a 1+2 reaction with only the elas- 
tic channel open, and with only relative s-waves involved, 
was considered. The applicability of the method is not 
limited to this particular case. In this work we shall con- 
sider processes involving t > relative angular momenta, 
and we shall derive the integral relations for the general 
case in which more than one channel is open. The only 
limitation is that we shall restrict ourselves to energies 
below the breakup threshold. Above it infinitely many 
adiabatic terms are in principle needed to describe the 
breakup channel, and although the same procedure could 
be used to describe it, we leave this particular case for a 
more careful investigation in a forthcoming work. 

A different aspect is the applicability of the method to 
describe 1+N reactions with A = 1 + N > 3. In this 
case, the main difficulty is to obtain the N+l wave func- 
tion in the internal region and the iV-body bound state 
function describing the asymptotic configuration. With 
this information, the integral relations apply exactly the 
same way as for the 1+2 case, but replacing the bound 
dimer wave function by the corresponding bound iV-body 
wave function. The extension of the adiabatic expansion 
to describe more than three-particles is possible. The de- 
pendence of the hyperangular part consists in (3JV — 4) 
hyperangles and, in the case of systems of identical parti- 
cles, the problem of constructing a ^4-body wave function 
with the proper statistic has to be faced. First applica- 
tions of the HA expansion to describe a four-body system 
already appeared [17] . In this work, however, we restrict 
the discussion to 1 + 2 reactions. 

In section II we describe the details of the formalism, 
first describing the adiabatic expansion in a multichan- 
nel reaction, and second showing how the corresponding 
/C-matrix (or equivalcntly the 5-matrix) can be obtained 
from the asymptotic wave function. In section III the 
integral relations for the same multichannel reaction are 
derived. They permit to extract the JC- (or S-) matrix 
requiring only knowledge of the internal part of the wave 
function. The results are shown in section IV. In section 
II V Al we consider a test case with only the elastic channel 
open. We investigate a three-body process which is fully 
equivalent to a two-body reaction, for which the phase 
shifts can be easily computed. This can then be used 
to test the accuracy of the integral relations method as 
well as the convergence pattern in the adiabatic expan- 
sion when I > partial waves are involved. In section 



IV B we investigate the elastic collision between a 4 He 
atom and the weakly bound ( 4 He)2 dimer. Finally, in IV 
C we apply the method to the collision involving a 6 Li 
and two 4 He atoms. In particular we shall consider in- 
cident energies such that the two possible incoming and 
outgoing channels, ( 4 He,( 4 He— 6 Li)) and ( 6 Li,( 4 He)2) are 
both open. The summary and the conclusions are given 
in section V. In appendix A we show the derivation of 
the Kohn Variational Principle for a multichannel pro- 
cess and, finally, in appendix B we have collected some 
technical details of the use of the integral relations when 
projected two-body potentials are employed. 



II. FORMALISM 

A. General features of the HA expansion 

In this work we consider a process where a particle hits 
a bound two-body system. We assume the incident en- 
ergy to be below the breakup threshold in three particles. 
This means that the total three-body energy E, which is 
the sum of the incident energy k 2 /2p > (p, being the re- 
duced mass between the incident particle and the dimer) 
and the two-body binding energy E2b, is negative. In this 
way only elastic, inelastic, and rearrangement processes 
are possible. 

The reaction under study is therefore a three-body pro- 
cess, which as usual, can be described through the x and 
y Jacobi coordinates: 



rrijUik 



m(m,j + TOfc 



(1) 



Vi = 



rrii (nij 



m k ) 



m(mi 



Vi - 



m k r k 



m k 



where and are the mass and coordinate of particle 
i and m is an arbitrary normalization mass. From the 
Jacobi coordinates one can construct the hyperspherical 
coordinates, which contain a radial one, the so-called hy- 
perradius p (p 2 = \J x 2 + yf) and the five hyperangles 
il ([il] = [ati,Cl x ,il y ]). The hyperangle on is defined as 
tan on = Xij yi and Q x and Q y give the directions of Xi 
and y i . The five hyperangles depend on the particular 
ordering of the particles chosen in the definition of the Ja- 
cobi variables. Three different sets are possible by cyclic 
permutations of the indexes i,j, k. In the following the 
Jacobi coordinates x and y and the corresponding hyper- 
angular coordinates are given using the natural ordering 
of the particles i, j, k = 1, 2, 3. 

Following Ref. Q we give a brief description of the HA 
method. In hyperspherical coordinates the Hamiltonian 
operator H takes the form: 
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where T p = J^p + is the hyperradial kinetic en- 
ergy operator, G 2 is the grand-angular operator and 
V(p, ft) — J2i Vi( x i) is the potential energy (i runs over 
the three Jacobi systems). 

The adiabatic expansion is based on the assumption 
that when describing a particular process, the hyperan- 
gles vary much faster than the hyperradius p. Under this 
assumption it is possible to solve the Schrodinger equa- 
tion (W — E)^> = in two steps. In the first one the 
angular part is solved for a set of fixed values of p. This 
amounts to solve the eigenvalue problem 
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for each p, which is treated as a parameter. 

The angular functions {$ n (p, Vl)} are used to construct 
the HA basis in which the basis elements form an or- 
thonormal basis for each value of p. The full three-body 
wave function is then expanded as: 
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J2f n (p)$ n ( P ,n). 
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Obviously the summation above has to be truncated, and 
only a finite number ha of adiabatic terms are included 
in the calculation. For simplicity we are omitting in \f r , 
f n , and $„ the quantum numbers giving the total three- 
body angular momentum and its projection. 

In a second step, the radial wave functions f n (p) in 
the expansion of Eq.Q are obtained after solving the 
following coupled set of radial equations: 
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where the operator T-L nn i acts on the radial functions and 
takes the form 
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for the diagonal terms, and 
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Hnn' = - ^ (^ P nn> (p) ^ + QnW (p)j (7) 

when n n' . 

The coupling terms P nn i and Q nn > in the expressions 
above follow from the dependence on p of the HA basis. 
Their explicit form is 

Pnn'(p) = ($„G>,n)|^|*n'G>,n)) n 



(p) = (* n (p,n)\^\$ n ,(p,n)) a , (8) 



where ()q represents integration over the five hyperangles 
only. 
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FIG. 1: Typical effective adiabatic potentials for a three- 
body system where two two-body bound states are present. 
The two lowest adiabatic potentials go asymptotically to the 



binding energies E^f and E 2b ' of the two-body bound states. 



For a given three-body energy E, when E 2b 

(2) 

one channel is open, while when E 2b 
are open. 



$<E< only 
< E < both channels 



The one-dimensional set of coupled differential equa- 
tions given in Eq.([S]) can be written in a matrix form 
as: 
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and the three-body wave function is 
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It is important to note that the diagonal terms T-L nn in 
Eq. ([6]) contain the angular eigenvalues A„ (p) introduced 
in Eq.(j3|). They appear in the effective adiabatic poten- 
tials, which are given by: 
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A typical behavior of the adiabatic potentials is shown 
in FigJTJ They correspond to a three-body system where 
two of the two-body subsystems have a bound state. This 
is reflected in the fact that the two lowest effective adi- 
abatic potentials go asymptotically to the binding ener- 
gies £^6 and E^ of each bound two-body system. The 
angular eigenfunctions associated to these two adiabatic 
potentials have the general asymptotic form [9J: 
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where for this particular case n = 1 , 2 and we have now 
made explicit the quantum numbers. The wave func- 
tion ip n *(x), normalized to 1 in the a;- Jacobi coordinate, 
describes the bound two-body system associated to the 

(n) 

effective potential V e > t , whose angular momentum is j x . 
Asymptotically it tends to the bound state wave function 
of the corresponding two-body subsystem. The spin func- 
tion Xs y describes the spin of the third particle, which 
couples to the orbital angular momentum l y (associated 
to the Jacobi coordinate y) to give total angular momen- 
tum j y . Finally, j x and j y couple to the total angular 
momentum J with projection M of the three-body sys- 
tem. 

The analytic form given in Eg. (1121) for the asymptotic 
expression of the angular eigenfunction $^ M (fi, p) makes 
evident that it describes an asymptotic spatial distribu- 
tion for the three particles corresponding to two of them 
forming a bound state, described by ip n x (x), and a free 
third particle moving in the continuum. In other words, 
the effective adiabatic potentials associated to angular 
eigenfunctions with the asymptotic form of Eq. (fT2")) are 
the ones describing the possible incoming and outgoing 
channels of a process where a particle hits a bound state 
formed by the other two. 

In Fig JT] the different regions defined by the energy of 
the incident particles are depicted. All the three-body 



energies E such that E% < E < E { 2b > (like E^ in the 
figure) correspond to processes where only one channel 
is open. Only the elastic collision between the third par- 
ticle and the bound two-body state with energy E^f is 
possible. When the three-body energy increases up to 
the region E^ < E < (E^ in the figure) a second 
channel is open. Two different collisions are now possi- 
ble, the one where a particle hits the bound state with 
binding energy E<u and the one where a particle hits 

(2) 

the state with binding energy E^b ^ n ^ ne same wa Y; 
each of these reactions has two possible outgoing chan- 
nels, corresponding to the two allowed bound two-body 
states and the third particle in the continuum. In par- 
ticular, in this energy range the rearrangement process is 
open. When E > the breakup channels are also open. 
They are described by the remaining infinitely many adi- 
abatic potentials. Processes with breakup channels open 
will be investigated in a forthcoming work. 

Therefore, for processes where no channels are open, 
the full three-body wave function has actually n differ- 
ent components. We shall denote them by corre- 
sponding to the process with incident channel i. Each 
of the three-body functions ^ is then expanded as in 
Eq.Q, but the radial functions need now an additional 
index i (i < no) indicating the incident channel to which 
they correspond: 
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For each open channel i (i = 1,2, •• • , no) the corre- 
sponding radial wave functions f n i(p) satisfy the set of 



radial equations of Eq. §Q , and Eq. © can be generalized 
to 
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which describes the full process. The full three-body 
wave function is now given by: 
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B. Asymptotics: /C-matrix and <S-matrix 

For scattering states and energies below the breakup 
threshold (E < 0), the Eqs. (fT4")) decouple asymptotically, 
and for a given incident channel i (i — 1, • • • , no) the only 
equations surviving are the ones of the form: 



(n nn -E)f m { P ) = o 



(n = ,n>A), (16) 



which, by use of Eq.(|6|), can be written as: 
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where V^l is given by Eq. fTTj) . 



When n corresponds to a closed channel, the radial 
wave functions f n i vanish asymptotically. When n cor- 
responds to an open channel, the asymptotic behavior of 
f n i is dictated by the asymptotics of the corresponding 
adiabatic potential V e fi . A careful analysis of the large 
distance behavior of the A„(p) and Q nn (p) functions in 
the case of bound two-body subsystems can be found in 
Ref. [9J]. In particular, Eqs.(91) and (93) of that reference 
allow to rewrite the above equation for the case n < no 
as: 
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where 



i?2 b is the binding energy of the bound two-body system 
associated to the open channel n, and £ y is the orbital an- 
gular momentum associated to the Jacobi coordinate y, 
which amounts to the relative orbital angular momentum 
between the projectile and the two-body bound target. 
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From Eq. ([18|) it is now clear that the asymptotic be- 
havior of the f n i functions (n, i < n ) is given by: 
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where ji and r]t are the usual regular and irregular 
spherical Bessel functions, respectively. The superscript 
(K) indicates that, with this particular choice, the co- 
efficients A-„ and will permit to extract the /C- 
matrix. Conversely, using the spherical Hankel functions 
in Eq. (|20|) the coefficients will form the 5-matrix and the 
superscript (S) will be used (see below). 

Therefore, asymptotically, the matrix containing the 
radial wave functions in Eq. ([15[) reduces to the no x hq 
matrix #)j+fiWy, where and B {K ^ are no x no 
matrices whose components are the A)*P and co- 



efficients of Eq. ([20]) . and J and Y are two no x no di 
agonal matrices with diagonal terms 
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respectively. Thus, the asymp- 



totic behavior of the full three-body wave function (fTS)) 
can be finally written as: 



where Fp and G p n> are column vectors with no 
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terms of the form (y ky n ^je y (k y n ^ p)^ n /p 3 ^ 2 
( V k y n) m y (ki n) p)$ n / P 3/2 ) , respectively. 



and 



From Eq. (f2lj) we then have that for a given incident 
channel i the asymptotic form of the corresponding three- 
body wave function ([T3")) takes the form: 



is the /C-matrix of the reaction, whose dimension is no x 
no (with no being the number of open channels) . 

The discussion in this subsection could have also been 
made by replacing jt v and ij£ y in Eq. ([20|) by the spherical 

Hankel functions \if^ and hf~\ respectively. This would 
then lead to: 
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where now F P (S) and G { P S) are column vectors with 
no terms of the form ( y k y n ^ h^f\k y n ^ p)^ n / p 3 ^ 2 ) and 



k y n) (k y n) p)<^ n / p 3 / 2 ). rcspccuvcly. We can then 
write: 



where 
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is the so called 5-matrix of the reaction. The S and 
JC matrices are related through the well known simple 
expression: 



S = {I + - iK)- 1 . 



(29) 



It is important to keep in mind that while A^ K \ B^ K \ 
and K, are real, the matrices A^ s \ B^ s \ and S are in 
general complex. 



III. SECOND ORDER INTEGRAL RELATIONS 
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and where we have made use of Eq. (fT2j) . which relates 
the angular eigenfunction and the two-body wave 

function ip^ . When two or three identical particles are 
present in the system, these functions should be correctly 
symmetrized or antisymmetrized depending on whether 
they are either bosons or fermions. 
From Eq. (|21[) we can now easily write: 



icg p k: 



where 
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In Ref . [14j the applicability of the HA expansion to ex- 
tract phase shifts for 1+2 reactions when only the elastic 
channel is open has been discussed. In that reference it 
was found that, when increasing the number of adiabatic 
channels ua included in the calculation as much as possi- 
ble, the difference between the computed phase shift and 
the exact value remains significant. As mentioned in the 
Introduction, this is related to the fact that the asymp- 
totic structure of the system has to be describe in terms 
of spherical Bessel functions depending on k y y, where y is 
the modulus of the Jacobi coordinate between the center 
of mass of the outgoing bound two-body system and the 
third particle. Instead, the asymptotic behavior using 
the HA basis is given in terms of spherical Bessel func- 
tions depending on k y p. Since the equivalence between 
k y y and k y p is not matched for any finite value of p, the 
correct boundary condition is only achieved at p rs oo 
and riA — > oo. 

For a general multichannel process the adiabatic ex- 
pansion obviously shows the same deficiency. The correct 
asymptotic wave function is given by Eq. (|21[) . but where 

F { P K) and G { P K) in Eq.® have to be replaced by 
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and G^ K \ which are column vectors whose n-th element 
is 
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In these expressions y n refers to the modulus of the 
Jacobi coordinate describing the center of mass of the 
bound two-body system ipf? and the third particle. 

It is important to recall that the Bessel functions r/g 
are irregular at the origin, which creates difficulties from 
the numerical point of view. It is then convenient to 
regularize such function, in such a way that G { n K) given 
in Eq.(|30p has to be replaced by: 
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where 7 is a non linear parameter. The results are stable 
for values of 7 within a small range around 7 ~ l/r , 
with ro the range of the potential. 

For simplicity in the notation, from now on, we shall 
refer to the matrices {F^ K \G^} as {F, G}, in such a 
way that we can write the asymptotic behavior of the 
wave function as: 



* -> AF + BG, 



(32) 



and K = -A~ X B. 

The vectors F and G satisfy the following normaliza 
tion condition: 
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where I is the identity matrix. In Eg. (1331) we have intro- 
duced a notation to be used from now on in which the 
overlap of two vectors is a matrix whose elements are, 
for example, ({F\H - E\G)) lJ = {F^U - E\G 3 ). The 
normalization condition allows to extract a first order es- 
timate of the matrices A and B from the scattering wave 
function as 
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Clearly, when ^ is an exact solution of (% — E)^> = 0, 
the above expressions reduce to the the following integral 
relations: 
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Explicitly, each matrix element Bij and Aij is given by 
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(38) 



which can be seen as the extension to multichannel scat- 
tering of the expressions valid for the single channel case. 
Now the same formula applies for each possible incoming 
channel described by ^ and each possible outgoing chan- 
nel whose asymptotic analytic form is given by a linear 
combination of Fj and Gj . 

As demonstrated in Refs. [3, EH for a single chan- 
nel process, the relation /C = —A~ 1 B computed using 
Eqs. (j37l) and ([38]) can be considered accurate up to sec- 
ond order when a trial wave function \l/ t is used. More- 
over the two integral relations of Eqs. (f3"Tj) and (f3"5|) can be 
directly derived from the Kohn Variational Principle. As 
shown in Appendix A, the matrix form of of KVP, neces- 
sary to describe a multichannel process, establishes that 
each matrix element of A~ X B 2 is a functional given by 
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which is stationary with respect to variations of the wave 
function. Taking into account the general asymptotic be- 
havior in Eq. (|32[) . we can write the full trial wave function 
schematically as: 



*t = * c + AF + BG, 



(40) 



with ^ c — > as p — > 00. Furthermore ^ c can be ex- 
panded in terms of a (square integrable) complete basis 
{bi,i = !,••■ , rn}: 



( * c < x ^ 


I 


#c,2 


1 






V *C«o j 


\ 



/ Cu C12 

C21 C22 

\ Cno 1 Cno 2 



Clm \ 
C2m 



/ 61 \ 

62 



\ b m ) 
(41) 

The variation of the functional with respect to the linear 
parameters Cij and with respect to the matrix elements 
of A~ X B leads to: 



(* c \H-E\V t ) = 
(G\H-E\* t ) = 



(42) 



When VP is replaced by ^>t, the second expression above 
and Ea.([3"5f result: 



9m 

A = - — (* t \H-E\G). 



(43) 



Replacing now Eq. (j40|) into (|39|) . and making use of the 
Eas.(|32"]). we also get: 

B 2 " d =B 13t + ^(F\H-E\* t ) T , (44) 

and, taking into account that B 1 is given by Eq. tjM]) we 
can then obtain the final result: 

2m 



B 2 " = J^(^t\H~E\F) 

/Til f 

A = -jjr(*t\H-E\G), 



(45) 
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FIG. 2: Scheme of the model reaction used to test the inte- 
gral relations. The light projectile 3 hits the dimer made by 
particles 1 and 2. Particle 1 is infinitely heavy and particles 
2 and 3 do not interact. 

which according to Eqs. ([2"5")) and (j2T))) permit to obtain 
the second order estimate of the /C-matrix or the 5- 
matrix, K 2 or S 2 , respectively. 

In practical cases, application of the integral relations 
given in Eq. (l4"5")l require the calculation of each individual 
matrix element Aij and -By, which relate each possible 
incoming channel i described by with each possible 
outgoing asymptotics given by Fj and Gj . Details about 
the calculation of these matrix elements are given in ap- 
pendix B, in particular for the case of two-body potentials 
projecting on the partial waves. 

The integral relations of Eq.([4"5"j) depend on the short 
range structure of the scattering wave function "Ft as F 
and G are asymptotically solutions of (H — E)F, G = 0. 
This property allows for different applications of the in- 
tegral relations, as discussed in Ref. [18|]. In the present 
work the interest is given in the study of the pattern 
of convergence of K. 2 in terms of the number of equa- 
tions tia considered in the description of 'Ft using the 
HA expansion. As we will see, increasing tia, both ma- 
trices A and B 2 slightly change, showing individually 
a very slow rate of convergence. Conversely, its rate 
JC 2 " = —A~ 1 B 2 ™ shows a pattern of convergence simi- 
lar to that one observed in a bound state calculation. 



IV. RESULTS 
A. Test of the method: A model 1+2 collision 

To test the method we have chosen a 1+2 reaction 
where the target dimer is made by an infinitely heavy 
particle and a light one, and where we consider a projec- 
tile interacting only with the heavy particle (see FigJ5]). 
In the collision particle 2 does not play any role, and the 
process is equivalent to a two-body reaction between par- 
ticles 3 and 1. Therefore the results obtained through the 
three-body calculation and the integral relations can be 
easily tested by means of a simple two-body calculation. 

In particular, we consider a two-body target made by 
two spin-zero bosons with masses 0.5m and 10 12 m (with 



TABLE I: Partial wave phase shifts 5e for different values of 
tia (number of adiabatic terms used in the expansion (0])). In 
the last row, the result using a two-body calculation is shown. 



n A 


5s 


s p 


s d 


1 


40.554 


0.6658 


0.0136 


2 


38.988 


0.6892 


0.0113 


3 


38.642 


0.6921 


0.0121 


5 


38.693 


0.6911 


0.0119 


8 


38.702 


0.6918 


0.0118 


10 


38.701 


0.6918 


0.0118 


two-body 


38.699 


0.6917 


0.0117 



m = 938.69461 MeV) interacting via a simple central 
potential given by 

V a (r) = -80e~ r2/1 ' 63 (46) 

where r is given in fm and the strength in MeV. This 
system has only one s-wave bound state with binding 
energy -6.2757 MeV. 

The projectile, which is chosen to have a mass of 
0.51m, does not interact with particle 2, while it does 
it with particle 1 through the gaussian potential: 

V 13 (r) = -30e- r2 ^ 6 \ (47) 

where again r is in fm and the strength in MeV. This 
potential is not able to bind particles 1 and 3. Finally, 
as described above, V23 = 0. 

We have chosen an incident energy of 3 MeV, which 
implies a total three-body energy of —3.2757 MeV. We 
are then below the threshold for breakup of the two-body 
target, and only the elastic channel is open. Therefore, 
B and A in Eq.(|45p are just numbers, and they are such 
that t&nSi = —B/A (note that the definition of B in here 
and in [13, (3 have opposite sign). 

We have computed the phase shift for this reaction for 
relative s, p, and d waves between the projectile and the 
target. The convergence of the expansion Q is shown in 
table [H where we show the phase shift for the different 
partial waves and for different values of tia , which is the 
number of adiabatic terms included in the calculation. 
As we can see, inclusion of 8 to 10 adiabatic potentials is 
enough to reach convergence for the three partial waves. 
Furthermore, the converged result agrees with the phase 
shift obtained from the two-body calculation describing 
the collision between particles 3 and 1. 

The efficiency of using the integral relations in Eq. (|45"|) 
is made evident in Fig. [3J where we show the partial wave 
phase shifts 81 as a function of ha- The solid line gives 
the results obtained from the integral relations (given 
in table [J) , and the thick dashed line shows the results 
extracted by direct comparison of the computed asymp- 
totic radial wave functions and the analytic expression in 
Eq. (f2T)]) . The thin dashed line indicates the phase shifts 
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FIG. 3: (color online) Phase-shift for s, p, and d partial waves 
for the model reaction in section IIV Al as a function of the 
number of adiabatic terms included in the calculation. The 
solid line is the result obtained through the integral relations, 
and the thick dashed curve has been obtained from the asymp- 
totic expression (|20p . The thin dashed line is the result ob- 
tained from a two-body calculation. 



obtained from a two-body calculation. As we can imme- 
diately see in the figure, the pattern of convergence of the 
phase shifts obtained from Eq. (j2"0l (thick dashed curves) 
is very slow. A simple extrapolation of these curves up 
to the correct value permits to foresee that the number 
of adiabatic terms needed to obtain accurate values of Si 
is far larger than the one needed when the integral re- 
lations are used. In fact, at the scale of the figure, the 
calculations with the integral relations are already for 
riA — 4 indistinguishable from the correct result (see also 
table U). 

As seen in table HI the d-wave phase shift is already 
rather small and therefore, at the considered energy, the 
cross section contributions from higher angular partial 
waves is negligible. In Fig. @]the differential cross section 
of the process with cumulative inclusion of one (solid), 
two (dashed), and three (dot-dashed) partial waves is 
shown. As can be seen, s-, p and rf-partial waves are 
enough to obtain a converged cross section of the pro- 
cess. In fact, the partial waves beyond the s-wave have 
a modest contribution, and the total cross section ap- 
proaches quite a lot the characteristic sin 9 function of 
the s- waves. 



s - wave 
s,p - waves 

s,p,d - waves 




1 2 
9 (rad) 



FIG. 4: (color online) Cumulative contributions to the differ- 
ential cross section as a function of the scattering angle 6 for 
the model collision in section [IV Al 



B. A realistic case: The 4 He - 4 He 2 collision 

In this section we discuss an interesting simple physical 
case, but technically similar to the model case described 
in the previous section. This is the collision of a 4 He atom 
into the weakly bound 4 He2 dimer. The helium dimer has 
a single s-wave bound state, and as soon as the incident 
energy is below the dimer breakup threshold, again only 
the elastic channel is open. 

The two-body helium-helium interaction is chosen to 
be the simple effective gaussian potential given in (l9j . 
This is enough for our purpose of illustrating how this 
kind of processes can be easily described by use of the 
integral relations. This potential is built to reproduce 
the s-wave scattering length (189.054 a.u.) and effective 
range (13.843 a.u.) of the LM2M2 interaction (20|, and 
it is given by: 



V 2B {r) 



-1.227 e- r2 / ia032 



(48) 



where r is given in a.u. and the strength is in K. This po- 
tential leads to a bound + 4 He2 dimer with a binding en- 
ergy E d = -1.2959 mK, a scattering length a = 189.947 
a.u., and an effective range of 13.846 a.u.. Simple repre- 
sentations of the atom-atom potentials are often used to 
describe reactions in the ultracold regime (see for exam- 
ple Refs. [H], HH). In this regime the process is largely 
independent of the shape of the potential and can be 
characterized only by the scattering length. 

With this interaction the helium trimer has two bound 
states at —150.0 mK and —2.467 mK. These states have 
been obtained using the gaussian potential active only in 
s-waves. Increasing the number of partial waves up to 
P-x = iy = 8 results in a very small change for the ground 
and excited state binding energies, which become now 
150.4 mK and —2.472 mK, respectively. In fact, more 
than 99% of the norm of the bound state wave functions 
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TABLE II: Partial wave phase shifts Si are given showing 
convergence of these values with the number ua of adiabatic 
potential used in the calculation. The last row shows the 
result when the Hyperspherical Harmonic method [J] is used. 



riA 


6 S 


dp 


s d 


Sf 


1 


-39.72 


-13.19 


2.01 


-0.27 


2 


-40.30 


-13.13 


2.11 


-0.28 


4 


-40.43 


-13.11 


2.13 


-0.28 


8 


-40.50 


-13.11 


2.14 


-0.28 


18 


-40.54 


-13.11 


2.14 


-0.28 


22 


-40.54 


-13.11 


2.14 


-0.28 


HH-calculation 


-40.55 









is provided by the lowest adiabatic term, whose corre- 
sponding adiabatic potential is close to identical in both 
calculations. Accordingly, in the following we restrict the 
calculations to include only the i x = channel. 

When the LM2M2 potential is used, these two states 
are found to have binding energies —126.4 mK and 
—2.265 mK, respectively |8j. As we can see, the ground 
state is not very well reproduced when the gaussian ver- 
sion of the potential is used. In this very deep state, the 
three atoms are close to each other and the correct struc- 
ture can not be described with the simplified potential. 
Conversely, the excited state which has the characteris- 
tic of an Efimov state has an structure in which the third 
atom orbits very far from the bound state of the other 
two. This particular structure is well described by the 
attractive gaussian potential. 

In order to study the convergence properties of the HA 
expansion for a 1 + 2 collision, we have chosen an incident 
energy of 0.5 mK (or a three-body energy E — —0.7959 
mK) . The phase shifts for the different partial waves have 
been computed as in the previous subsection. The results 
are shown in table |TT] for s, p, d, and / waves. A good 
convergence is obtained already after inclusion of about 
10 adiabatic terms, except for s-waves, where about 18 
are needed. The last row in the table shows the phase 
shift obtained for an s-wave collision when the Hyper- 
spherical Harmonic method is used. The two methods 
are in close agreement. 

At this particular energy we have calculated the differ- 
ential cross section. Fig. [5] shows the cumulative contri- 
butions of the s, p, d and / partial waves. We observe 
that the p-wave contribution is rather important and pro- 
duces a deviation from the sin# shape. Moreover, four 
partial waves are needed to reach a good convergence. 

In Fig. [6] we show the computed s-wave phase shift as a 
function of the incident energy (E — Ed). Our results are 
given by the stars. For comparison we also show the re- 
sults reported in [23[, (24j, and [25[ (solid curve, squares, 
and dashed curve, respectively). As we can see, the phase 
shifts obtained in this work are a few degrees above the 
ones obtained in the previous calculations, where the 4 He 
- 4 He interaction is treated more in detail. In fact the rea- 
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FIG. 5: (color online) The same as in Fig|4]for the 4 He - 4 He2 
collision. With four partial waves the convergence is fairly 
good. 
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FIG. 6: (color online)s-wave phase shift 5 a as a function of 
the incident energy for the 4 He - 4 He2 collision. The stars are 
the results obtained in the present work. The solid curve, the 
squares, and dashed curve correspond to the results given in 
[2g|, 0! and respectively. 



son for this discrepancy is the hard core rep ulsion present 
in the 4 He - 4 He interactions used in |23l - [25| . For the same 
reason the atom-dimer scattering length a a -d obtained 
with the gaussian potential used in this work, a a -d = 166 
a.u., differs from the typical values of around a a -d = 220 
a.u. (w 116 A) obtained when hard core potentials are 
used tHIH. 

C. A multichannel collision: The 4 He- 4 He- 6 Li 
system 

In this subsection a reaction where more than one 
channel is open is discussed. To this aim, we have chosen 
a process involving two helium and one lithium atoms. 
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The cross section for this kind of reactions is the nec- 
essary ingredient to obtain the recombination rate for 
such three-body systems. As quoted in [271 ] ■ where the 
three-body recombination for cold helium-helium-alkali- 
metal systems is investigated, such collision processes are 
important in ultracold gas experiments using buffer-gas 
cooling, since it might limit the lifetimes of the trapped 
atoms. 

To describe this three-body system we take the same 
helium-helium interaction as in the previous section, 
which leads to a + 4 He2 dimer with a binding energy 
of —1.2959 mK. The lithium-helium interaction is also 
chosen to have a gaussian shape, and it is taken to be: 

^Li-^Hc) M = -0-27368 e - 2 /20.i4 2 ^ (4Q) 

where r is in a.u. and the strength is in K. The pa- 
rameters have been adjusted to give a scattering length 
of —173.5 a.u. and an effective range of 26.475 a.u. in 
agreement with the values obtained in [28[ (a = —173.8 
a.u. and r e =26.483 a.u.), where the more sophisticated 
KTTY potential is used. This potential leads to a + 
bound 6 Li- 4 He system with a binding energy of —1.4225 
mK. 

The adiabatic potentials obtained for the 4 He- 4 He- 
6 Li three-body system follow the same pattern as the 
potentials in Fig[TJ where E^ corresponds now to the 
binding energy of the 6 Li- 4 He dimer (—1.4225 mK) and 
corresponds to the binding energy of the 4 He2 
dimer (—1.2959 mK). The three-body system presents 
one bound state at E = —58.12 mK. 

Thus, as soon as the three-body energy lies in the same 
region as E^ in the figure, two different channels are 
open. One of them corresponds to a bound 6 Li- 4 He dimer 
and the second 4 He atom in the continuum (we shall refer 
to it as channel 1), and the other one corresponds to the 
bound 4 He2 dimer and the 6 Li atom in the continuum 
(we shall refer to it as channel 2). In other words, when 
taking channel 1 as the incoming channel we are consid- 
ering a process were the 4 He atom hits a bound 6 Li- 4 He 
dimer, while when choosing channel 2 as the incoming 
channel we are then considering the process of a 6 Li atom 
hitting a 4 He2 dimer. For each of the two possible incom- 
ing channels we have two different outgoing channels, the 
elastic one and a rearrangement process where the pro- 
jectile is captured by one of the constituents of the dimer, 
while the second dimer constituent is released. 

The existence of two open channels implies that the 
JC- matrix (or the 5-matrix) is a 2x2 matrix, that can 
be obtained through the 2x2 matrices A and B 2 " d in 
Eq. (|45| . Each of the four terms in A and B 2 can be 
obtained as in Eqs. (f3"5|) and (|3"T|) where is the trial 
three-body wave function for the incoming channel i. Fj 
and Gj are the asymptotic functions given in Eqs. (f3"0)) 
for the outgoing channel j. These functions have to be 
symmetrized when the outgoing channel is 1, since the 
helium atom in the dimer is identical to the one moving 
in the continuum. For outgoing channel 2 this is not nec- 



TABLE III: /C-matrix elements are given as a function of the 
number of adiabatic potentials used in the calculation (ha)- 



riA 






i- ■ 

A-21 


/V22 


o 
I 


O Add 


-U.DOU 


-0.048 


1 /ill 
-1.411 


3 


-2.765 


-0.821 


-0.801 


-1.496 


4 


-2.691 


-0.775 


-0.776 


-1.468 


6 


-2.699 


-0.781 


-0.781 


-1.471 


8 


-2.702 


-0.783 


-0.783 


-1.471 


10 


-2.710 


-0.787 


-0.787 


-1.473 


14 


-2.714 


-0.790 


-0.789 


-1.474 


18 


-2.712 


-0.791 


-0.790 


-1.474 



essary, since the dimer wave function in (1301) is already 
properly symmetrized. In practice, the symmetrization 
gives out a factor of y/2 in Eqs. (l3T|) and (|3"51) when j=l. 
The calculation of each of these terms is formally identi- 
cal to the case with only the elastic channel open. 

In the calculation here we have chosen a three-body en- 
ergy of —0.7959 mK, which represents an incident energy 
of 0.6266 mK when channel 1 is the incoming channel, 
and 0.5 mK when channel 2 is the incoming channel. For 
simplicity we restrict in this section to relative s-waves 
between the projectile and the dimer target. The com- 
puted result for the four terms of the /C-matrix are shown 
in table IIHI for different values of ua ■ As seen in the ta- 
ble, again a reduced amount of adiabatic terms permits 
to reach a reasonable convergence in the /C-matrix. 

From the computed /C-matrix we can now easily obtain 
the iS-matrix by means of Eq. ([29| . This leads to S\\ = 
-0.673 - 0.663i, S 12 = S 21 = -0.285 4- 0.162i, and S 22 = 
—0.224 — 0.918i The square of these elements, |>%| 2 , 
indicates the probability for the process with incoming 
channel i to end up in channel j. In this particular case 
we get \Sn | 2 = \S 22 1 2 = 0.892 and |Si 2 | 2 = |S 2 i| 2 = 
0.108. 

It is important to note that the matrices A and B 2 
are not unique. A different definition of the normaliza- 
tion of the asymptotic states would result into new A 
and B 2 " matrices which would obviously lead to the 
same /C-matrix. In particular, A and B 2 " do not fulfill 
the property of being symmetric, but they lead to a K- 
matrix with the correct hermitian condition. Moreover, 
using Eq. (|29|) , the computed 5-matrix automatically sat- 
isfies the unitarity condition S^S = I. 



V. SUMMARY AND CONCLUSIONS 

In this work we have discussed the general form of 
the integral relations that were introduced in [l4|, fl8| . 
These relations are derived from the Kohn Variational 
Principle and they permit to exploit the particularities 
of the adiabatic expansion method to describe scattering 
states. In particular, in [l4|, [U| it was shown that the 
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convergence of the computed scattering phase shifts in 
terms of the adiabatic terms included in the calculation is 
rather fast. The convergence pattern results to be similar 
to the one of a bound state calculation. The reason for 
this success is that when using the integral relations only 
the internal part of the wave function is needed, and an 
accurate calculation of it requires a smaller amount of 
adiabatic terms than when computing the wave function 
in the asymptotic region. 

The applications given in 0, HH were limited to pro- 
cesses involving only relative s-waves and with only one 
channel open. In this work we have explicitly derived the 
integral relations from the KVP in the case of multichan- 
nel reactions and we have computed phase-shifts up to 
/-waves. Furthermore, we have used a vectorial notation 
for the wave function such that all the possible channels 
are simultaneously represented. With this notation the 
coefficients weighting the regular and irregular part of 
the asymptotic wave function are uq x uq matrices (with 
no being the number of open channels) and each term of 
these two matrices is obtained from an integral relation. 
Finally, the /C-matrix for a given process is obtained as 
the product of two no x tiq matrices. 

Although the method derived is completely general, in 
this work we have restricted ourselves to describe 1 + 2 re- 
actions with projectile energy below the breakup thresh- 
old in three outgoing particles. Therefore, only elastic, 
inelastic, and rearrangement processes are possible. 

To test the method when including relatives partial 
waves higher than zero, we have first used a toy model 
such that the three-body reaction is fully equivalent to 
a two-body process. In this way the correct phase shift 
can be easily computed through a simple two-body cal- 
culation. We have found a slow convergence of the phase 
shifts when extracted from the asymptotic part of the 
radial wave functions. Conversely, the rate A~ X B con- 
verges much faster and the result stabilize with a rather 
small number of adiabatic channels. The convergence is 
equally fast for all the partial waves, and around 10 adia- 
batic terms are enough to reach a good convergence. Fur- 
thermore, the phase shifts obtained with the two-body 
calculations are well reproduced. 

As the next step, we have analyzed a more physical 
case, in particular the He - 4 He2 collision. Since we have 
considered energies below the 4 He2 breakup threshold, in 
this reaction only the elastic channel is open. This is a 
process technically analogous to the previous schematic 
case for which the method has been proved to work. In 
fact, a similar pattern of convergence is found for the 
different partial waves included in the calculation. Inclu- 
sion of partial waves with £ up to 3 are needed to obtain 
a converged cross section for the process. For s-waves the 
computed phase shift reproduces the one obtained with 
the Hyperspherical Harmonic expansion method. 

Finally, we have considered a process with two open 
channels. We have chosen a three-body system made by 
two He and one 6 Li atoms, where two different dimers, 
4 He2 and 4 He - 6 Li, are possible. We have therefore simul- 



taneously investigated the collision between a 4 He atom 
and a 4 He - 6 Li dimer, and the one between a 6 Li atom 
and a 4 He2 dimer. For both reactions two possible outgo- 
ing channels (elastic and rearrangement) are permitted. 
We have then used the method to obtain the 2 x 2 JC- 
matrix. Again, we have found a fast convergence of the 
four terms in K,. Furthermore, the computed /C-matrix 
satisfies the required hermitian condition, as well as the 
fact of leading to a unitary 5-matrix. 

Summarizing, we have shown that the integral rela- 
tions can be easily applied to reactions involving non- 
zero partial waves and with more then one channel open. 
Also in this case, the hyperspherical adiabatic expansion 
method is a highly efficient tool that permits to obtain 
scattering wave functions. The /C-matrix, and therefore 
also the 5-matrix, converges rather fast. Also, since the 
hyperspherical adiabatic expansion method permits to 
identify every single incoming and outgoing open chan- 
nel with a single adiabatic term, the dimension of the 
matrices to be computed is rather modest, typically of 
the same size as the number of open channels in the re- 
action. 



Appendix A: Matrix form of the Kohn Variational 
Principle 

This derivation is completely analogous to the deriva- 
tion presented in Ref. [29|]. The only difference is that, 
in order to represent each possible incoming channel, we 
use the vectorial notation for the wave functions as intro- 
duced in the present work, where the total wave function 
\& has the form given in Eq. (fT5|) . 

We start by taking the matrix given by: 

I=(^f\U-E\9) (Al) 

which vanishes when ^ is the exact wave function. We 
then introduce a test wave function *f? t — ^ + so that 
its radial wave functions verify: 

/™(0)=0 (A2) 

f m (p -+ oo) -J- Jkf ] p (A m]ly {k^p) + S^(fcWp)) ■ 

We can then write = f n i + fifni, where Sf„i satisfies 
that 

Sfm(p) -> 4W ] R m v (ki n) P ) SB m . (A3) 

The matrix in Eq. (|Al[) . evaluated at the test wave 
function, is: 

I t =8I={9t\H-E\^t) (A4) 

Using now that the exact wave function verifies that 
(H — E)\^) = 0, and keeping only the first order terms, 
the matrix above can be written as: 

51= (V\H - E\6V) - (Sty\H-E\y) T (A5) 
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Using the expansion of \E' given in Eq. (IT3|) and the 
analytical expression of the operator H from Eq. ((2J, 
it can be seen that each matrix element of 61 takes the 
form: 



SL, 



n 

m ^ 



' d (df ni d(Sf nj ) 



dp-riPiMfmSfm), (A6) 
o a P 



where Py (p) are the coupling terms appearing in Eq. (J7J , 
and given in Eq.®, which vanish when p tends to zero 
or to infinity. Therefore the last term in the previous 
expression vanishes, and, since /m(0) = 0, we get: 



61, 



2m 



E 



dfn 



■6f, 



d(6f n 



' fn 



dp dp 
which, using Eq.® and Eq.(|A"3")). leads to 



SL. 



2m 



Y, A ™SBj 



2m 



or, in a more compact way: 

6[I+^AB T 
2m 



(A7) 



(A8) 



(A9) 



Since for the exact wave function ^> we have that 1 = 0, 
we finally get: 



2m 



-AB 



It 



2m 



-AB 



(A10) 



which becomes a variational principle for AB T . There- 
fore, given a test wave function \l/t, we obtain a second 
order correction for AB T as: 



A(B ) T = AB 1 



2m 



9 t \H-E\9t 



(All) 



If we now multiply from the left by A -1 and from the 
right by (^4 _1 ) T , and make use of the fact that the K,- 
matrix is symmetric, i.e. A~ 1 B = _B T (^4 _1 ) T , we then 
finally get the expression given in Eq. (1551) . 



Appendix B: Calculation of the integrals in A and B 

In this appendix we give details of the calculation of 
the integrals A and B, in particular when using two-body 
potentials projecting on partial waves. To this aim let us 
start from the general expression for A and B in Eq. P5l) , 
and write the C operator in its explicit form 



2m , - . 2m / h 2 9 



(Bl) 



2m 



Vxi +Vi(x 1 ) + V 2 (x 2 ) + V 3 {x 3 ) -E d -Ek 



where Ed is the binding energy of the dimer and Eq is the 
incident energy of the projectile. The Jacobi coordinates 
{x\iUi) are defined such that x\ connects the two par- 
ticles in the dimer (particles 2 and 3). The coordinates 
x 2 and x 3 are related to the distances between particles 
1 and 3, and between particles 1 and 2, respectively. 

Using F and G as defined in (|30|) we can rewrite 
Eq.® as: 



B 2 " 
A 

where 



Z777 * 

- ¥ -<*\V 2 (x 2 ) + V 3 (x 3 )\F(x 1 ,y 1 )> (B2) 

2m ~ ~ 

-— < *\V 2 (x 2 ) + V 3 {x 3 )\G{x x , yi ) > +7 V 



h = - < *l V«i -k 2 m \G{x u yx) >, 



(B3) 



and where G refers to the regularized function (|3ip . 

If we call B 2 = Ib and A = I a + -fy> we nave that, 
after substitution of Eq. (|3U| , the integrals in (|B2[) can be 
written as well as: 

2m f 

I am = E / dpp^Q^ip^Viixi) (B4) 



h 2 ^ 

i=2,3 



9aj3{Vi, kyi) (asi) ® F V (^i)] 



where 



k y {kyV) 
-Ve v {k v y) (1 



(B5) 



7?A f y + 1 



and where for simplicity in the notation we have assumed 
that the particles have zero spin. The corresponding ex- 
pressions in this appendix for particles with spin will 
follow immediately by coupling the orbital part in the 
expressions above to the corresponding spin part. 

If the potential operator is given as a sum of projectors 
on partial waves, we have that: 



E 



VL (xi)\e xi m e x .){£ Xi m e 



(B6) 



where V^ x represents the interaction between particles 
j and k when they are in a relative partial wave with 
angular momentum £ Xi . 

If we also consider Eq. (UJ) and expand the angular func- 
tions $ n (P) in terms of the hyperspherical harmonics 

($„(p,n)^ = EKi ? i y CP w ( P )yfJ y M (m we can 
then obtain the following expression for the potential op- 
erator acting over the three-body wave function: 



Ki 



Expanding now the hyperspherical harmonics in terms 
of the Jacobi polynomials /A+VVv+i/a [K = 2v + £ x + 



13 



£ y ) with normalization coefficients N^ ly (sec [9] for de- 
tails) , we have that the the integrals I a and 1b can then 
be explicitly written as: 



Ia,b = E I dpp 5 (sin ai) 2 (cos arfdat-^ 



(B8) 



E/«W E E 41^4^ (Binoi^Coosoi)^ 
^ +1/2 ^ +1/2 (2cos a 0(^ M |^(^0l« M 



wave with angular momentum £ x . In general, for par- 
ticles with spin, the partial waves are identified by the 
quantum numbers {£ x , s x , j x }, where s x is the coupling 
of the spins of the two particles, which in turn couples to 
i x to give the total two-body angular momentum j x . In 

this case the matrix element lyf ~ M \Vt (x)\y¥^ M ^ 

Eqs. p7| or (|B8|) has to be replaced by: 



^/ 1 ( a; i)5A,fl( fc »n2/i) 



where is the radial part of the dimer wave function 

It is important to note that in Eqs. (|B2p and (|B4[) the 
potential operators and the functions F and G are writ- 
ten in a different Jacobi set. Therefore, when computing 
the integrals one has to rotate the whole integrand into 
the same Jacobi set. This is made in the expression above 
by the function Rij, which is a rotation function defined 
as: 



V J 3 w £xj£vj ( Xj , yj ) 



(B9) 



LM L 



LAIt. 



which rotates any function Wi x .i y (xj,yj) written in 
terms of the coordinates and angular momenta defined 
in the Jacobi set j into the coordinates and angular mo- 
menta corresponding to the Jacobi set i. 

As already mentioned, Ve(x) is the total two-body in- 
teraction when the two particles are in a relative partial 



y 



KL 



JM 



\Ve B s x jA x )\ 



y 



KL 



JM 



where s y is the spin of the third particle, \i 



(BIO) 
is three- 



body spin wave function and J, M are the total three- 
body angular momentum and its projection. In general, 
the partial wave two-body potential Vt xSm j m ( x ) could con- 
sist in a sum of central, spin-orbit, spin-spin and tensor 
potentials. Therefore, in this case, calculation of the ma- 
trix element in (|B10j) implies calculation of the matrix 
element of the corresponding spin-spin, spin-orbit, and 
tensor operators. In the simplest case with only a central 
potential and particles with zero spin the matrix elements 
of the potential operator reduce to: 

(y£g'\vt.w\y£z<) = yMs e ^j u j K , K , 

(Bll) 

simplifying the expression (|B8|) . 
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